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ABSTRACT 

We compute the corrections to the cosmological hydrogen recombination history due to delayed feedback of Lyman-series photons 
and the escape in the Lyman-continuum. The former process is expected to slightly delay recombination, while the latter should allow 
the medium to recombine a bit faster. It is shown that the subsequent feedback of released Lyman-n photons on the lower lying 
Lyman-(n - 1) transitions yields a maximal correction of AN e /N c ~ 0.22% at z ~ 1050. Including only Lyman-/? feedback onto the 
Lyman-tt transition, accounts for most of the effect. We find corrections to the cosmic microwave background TT and EE power 
spectra with typical peak to peak amplitude \ACf T /Cf T \ ~ 0.15% and \ACf E /Cf E \ ~ 0.36% at / < 3000. The escape in the Lyman- 
continuum and feedback of Lyman-o- photons on the photoionization rate of the second shell lead to modifications of the ionization 
history which are very small (less than |A_/V C /A/ C | ~ few x 10~ 6 ). 
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1. Introduction 

The tremendous advances in observations of the cosmic mi- 
crowave background ( CMB) temperature and polarization 
angular fluctuations (e.g. Page et al. 2006; Hinshaw et al.H2006h 
and the prospects with the Planck Surveyof] have motivated 
several groups to re-examine the problem of cosmological 
hydrogen recombination, with the aim of identifying previ- 
ously neglected physical processes which could affect the 
ionization history of the Universe at the level of > 0.1%, 
and may lead to modifi cations of the cosmological hydrogen 
recombination spectrum |Dubrovichll975l:lBeigman & Sunvaevl 
1 1978 1 : iRvbicki & dell' Antonio! 119931; IDubrovich & Stolvaroyl 
19951; iBurginl 120031; IDubrovich & Shakhvorostoval 12 004 ; 
Kholupenko et al. 2005; Wong et al. 2006; Rubi no-Martin et al l 



20061 IChluba & Sunvae v 2006a; IChluba et all 120071) . which 
could become observable in the future. For example, effects 
connected with the two-photon transitions from high s and 
d-states to the gro u nd state (|D~ubrovich & Grachev! [2005; 
IWong& Scottl 120071; IChluba & Sunvaevl 120071) the induced 
2s-two-photon decay ((Chluba & Sunvaev 2006b), the increase 
in the 2s-two-photo n absorption rate due to the la rge Lyman-o' 
spectral distortion dKholupenko & Ivan chik 2006), and details 
in the evo lution of the populations of the angular mom entum 
sub-states dRubifio-Martm et al.ll2006l IChluba et ai1l2007l) were 
discussed. All these studies show that sub-percent-level correc- 
tions to the ionization history do exist, which in pri nciple could 
bias t he values of the key cosmological parameters dLewis et al.l 
2006). It is clear, that when reaching percent-level accuracy in 
the determination of the key cosmological parameters or when 
considering signatures from inflation, e.g. the possibility of a 
running spectral index, accurate theoretical predictions of the 
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ionization history are required. Here we examine the effects due 
to feedback of hydrogen Lyman-series photons and the escape 
of photons in the Lyman-continuum on the ionization history 
during the epoch of hydrogen recombination, and the impact of 
this process on the Lyman-a distortion of the CMB blackbody 
spectrum. 

The strongest distortions of the CMB spectrum arising dur- 
ing the epoch of hydrogen recombination are due to the Lyman-q r 
transition and the 2s two-photon decay dZeldovich et al.| [l968; 
iPeeblesI 11968). The feedback of these excess photons on the 
phot oionization rates of the second shell has been considered 
by ISeager et al.l dl999i [2000), with no significant changes for 
the number density of free electrons, /V e . We confirm this state- 
ment for the feedback from Lyman-o' photons and find a max- 
imal correction of AN e /N e ~ 10~ 6 at z ~ 500, the redshift 
which roughly corresponds to the time, when the maximum of 
the CMB spectral distortion due to the Lyman-a transition is 
reaching the Balmer-continuum frequency. 

Ho wever, it has been shown by Kholupenko & Ivanchik] 
(2006) that the huge excess of photons in the Wien-tail of the 
CMB due to the Lyman-a distortion leads to an increase of 
the Is — > 2s two-photon absorption rate, which delays recom- 
bination and introduces corrections to the ionization history at 
the percent-level. Similarly, one expects some feedback of es- 
caping Lyman-series photons on lower lying Lyman-transitions 
and eventually the Lym an-qr transition. Recentl y, this process 
was also considered by ISwitzer & Hiratal d2007l) for both he- 
lium and hydrogen recombination, using an iterative approach. 
Here we directly compute the escaping Lyman-» radiation within 
the Sobolev approx imation (e.g. see iRubino-Ma rtm et al.| [2006; 
Chlu ba et al.l 12007) and include the feedback on the Lyman- 
(n - 1) transition at all times by evaluating the distorted spec- 
trum at the corresponding frequency. In addition, we assume that 
due to the huge optical depth in all the Lyman-lines, at the end 
only the Lyman-Q- line remains and all the Lyman-n lines will 
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Table 1. Appearance of the first few Lyman lines for the com- 
putations not including any feedback. In particular we give the 
approximate redshift, z max , at which AI v (z em ) is maximal (see 
Fig. [TJ|, and the redshift at which the peak of the line reaches 
the next lower Lyman resonance, zt ~ Zmax VLy(«-i)/vLy«, where 
VLyjt is the resonance frequency of the Lyman-A: transition. Also 
we give the approximate total number of escaping photons per 
hydrogen nucleus for each transition. 



Line 


n 


VLyn [HZ] 




Zf 


N y /Nn 


Lyman-ff 


2 


2.4674 x 10 13 


1400 




4.3 x 10- 1 


Lyman-/? 


3 


2.9243 x 10 15 


1450 


1223 


1.8 x 10 -3 


Lyman-y 


4 


3.0842 x 10 15 


1461 


1385 


2.9 x 10~ 4 


Lyman-(5 


5 


3.1583 x 10 15 


1464 


1430 


1.3 x 10~ 4 


Lyman- e 


6 


3.1985 x 10 15 


1467 


1449 


7.9 x 10~ 5 



be completely re-processed by the closest lower lying Lyman- 
(n - 1) transition. Therefore, feedback always only works as 
Lyn — > Ly(n - 1). It is also important to ensure that the line 
is not producing any feedback on itself. 

In addition, we estimate the escape probability for the 
Lyman-continuum and include the approximate net rate for di- 
rect recombinations into our multi-le vel code. As argued earlier 
(IZeldovich et al.lfT968tlPeeblesll968h . we find that the modifica- 
tion of the ionization history due to this process is indeed com- 
pletely negligible. 

2. Re-absorption of Lyman-series photons 

Although in computations with no feedback, the contributions 
of the np — » Is transitions from shells wi th n > 2 to the total 
Lyma n-series are rather small (see Fig. 8 in Rubino-Mart m et all 
2006), in the Wien-tail of the CMB the intensity of these lines is 
significantly larger than the blackbody (see Fig. [TJ. Therefore 
one expects that after redshifting each Lyman-transition from 
n > 2 will increase the absorption rate for the closest lower 
Lyman-series transition (n' = n — 1 > 2) and thereby slightly 
delay hydrogen recombination. For example, the Lyman-/? tran- 
sition peaks at z em ~ 1450 with contribution AI V /B V ~ 64% 
relative to the CMB blackbody spectrum B v . These excess pho- 
tons reach the central frequency of the Lyman-a transition at 
Z ~ 1223 and increase the Lyman-or absorption rate, thereby 
delaying recombination (see Table [TJ. It is easy to estimate the 
total number of photons released in the Lyman-/? transition us- 
ing Fig. [TJ Its value is close to ~ 0.18% of the total number of 
hydrogen nuclei (see Table [TJ for the other transitions). At the 
time of maximal feedback (zt ~ 1223) the number of free elec- 
trons is roughly 38% of the total number of hydrogen nuclei. 
Therefore the maximal effect one can expect for the feedback- 
induced corrections to the ionization history from the Lyman-/? 
line is ANJN e < 0.5%. 

In Fig. [2] we present the results for detailed computa- 
tions of feedback-induced changes in the free electron fraction. 
Including only Lyman-/? feedback on the Lyman-or transition ac- 
counts for most of the effect, leading to a maximal difference of 
ANJN e ~ 0.25% at z ~ 1040. 

When including in addition the feedback of Lyman-y on 
Lyman-/?, the maximal relative difference decreases to AN e /N e ~ 
0.20%, now peaking at z ~ 1070. This behavior is expected, 
since the Lyman-y photon reduces the transition rate in the 
Lyman-/? channel such that its feedback on the Lyman-or line 
should become slightly smaller. In addition, the number of elec- 
trons reaching the ground state via the Lyman-/? transition will 



io- 2 V 




-29 _ / ' ' ' 

10 e / / / Lyman-a v | 

: / // Lyman-P : 

.30 / / // Lyman-y 

^ z / / // Lyman-5 = 

: / //' Lyman-£ : 

,„-3i t l iz lA i i i i i i i i i 

600 700 800 900 1000 1100 1200 1300 1400 1500 1600 1700 1800 1900 




v [GHz] 



Fig. 1. Spectral distortions due to the first few Lyman-transitions 
at z — 0. No feedback Lyn — » Ly(n - 1) has been included. 
Upper panel: as a funct ion of the redshift of emission (see Eq. 1 
in Rubino-Martm et al. 2006). Lower panel: relative to the CMB 
blackbody as a function of observing frequency v. 



be reduced, which very likely leads to the small increase of the 
correction at z > 1200. 

Subsequent inclusion of more feedback results in an alter- 
nating behavior of AN e /N e . We found that for 30-shells, conver- 
gence is reached when including feedback from np-levels with 
n < 15. In Fig. [2] we also give the full result including all feed- 
back within a 30-shell hydrogen atom. In addition we found that 
the correction barely depends on the total number of included 
shells. Therefore we expect that even for computations with up 
to 100-shells the total feedback-induced relative difference will 
not change significantly. 

Examining the final differences in the Lyman-a distortion 
and the feedback-induced modifications of the Lyman-/? and y 
lines before their absorption within the corresponding lower ly- 
ing resonance (see Fig. [3]l shows that the situation is a bit more 
complicated. Including only Lyman-/? feedback the amount of 
emission in the Lyman-ff line reduces by ~ 0.8% at z ~ 1232. 
Note that this redshift is very close to the expected value, zt ~ 
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Fig. 2. Changes in the free electron fraction due to inclusion of 
Lyman-series feedback relative to the reference model without 
feedback. The computations were performed for a 15-shell atom 
where for the Lyman-series the escape of photons was modeled 
using the Sobolev approximation. The curves are labeled accord- 
ing to the sequence of feedback that was included. For example, 
4 — > 3 — > 2 indicates that Lyy — > Ly/3 and Ly/3 — » Lya feed- 
back was taken into account. Also the full result for the 30-shell 
hydrogen atom is shown, where feedback up to n = 30 was in- 
cluded. 



1223, for the maximal feedback (see TableQ]). On the other hand, 
due to the feedback-induced small changes in the populations 
of the levels, the net 2s- Is two-photon transition rate also in- 
creases with a maximum of roughly ~ 0.5% at z ~ 1034 (see 
Fig. ID, which partially cancels the delaying feedback effect on 
the Lyman-a transition. 

When including Lyman-y feedback, as expected, the emis- 
sion in the Lyman-/? line reduces, with a maximal difference at 
the expected value z ~ 1384. Adding Lyman-<5 feedback reduces 
the strength of the feedback of Lyman-y on the Lyman-/? transi- 
tion by a significant fraction. Comparing the relative strength of 
the Lyman-5 at its maximum with the Lyman-y line at zi (Fig.[TJ 
yields A/*(z max )/A/y (zi) ~ 1 /2, which indicates that Lyman-(5 
should be able to affect Lyman-y strongly. Looking at the rel- 
ative difference in the Lyman-y when including Lyman-<5 feed- 
back shows that at zi ~ 1430 the line is reduced by ~ 50%. 



3. The Lyman-continuum escape 

To estimate the maximal effect arising due to successful es- 
cape of photons from the Lyman-continuum we shall neglect 
any re-absorption of photons below the threshold frequency 
v c m 3.288 x 10 15 Hz. Although one can expect some modifi- 
cation due to the forest of lines in the Lyman-series close to the 
continuurrfl this approach will yield an upper limit on the effect, 
which turns out to be very small in any case. 

With this in mind, in the vicinity of the Lyman-continuum 
the photon number only changes by recombination and ioniza- 



2 Within a small energy-distance £ below the ionization potential 
these in principle can be treated as a c ontinuation of the pho toionization 
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tion to the ground state (we neglect electron scattering) 
1 <9M 

— A/ A/ III w r , 

r.2 



dt 



= N e N p f(T e ) cnJ^-e'^i - N u o- lsc N v . (1) 

Is C 

Here N e and N p are the free electron and proton number densi- 
ties, A^is is the number density of hydrogen atoms in the ground 
state and N v = I v /hv, where I v is the photon intensity. In addi- 

/ \LTE 

tion, f(T e ) = [N\ s /N e Np) and <x lsc is the ls-photoionization 
cross-section. As can be seen from Eq. ([T), photons are emitted 
with a spectrurrQ <p™ oc ct"i sc (v) y 2 e~ /,v / iB7 V . We now compute the 
probability that one single recombination will lead to a success- 
ful escape of photons. 

The optical depth for absorption of a photon, which has been 
emitted at redshift z em and is observed at frequency v ODS at red- 
shift z 0D s, within the Lyman-continuum is given by 



7c(Zem> Zobs? Vobs) 



Nu cr lsc (v z ) 



dz 



(2) 



H(l+z) 

where v z = v t, s (l +z)/(l + z bs)- If we substitute v = v z and use 
dv = v obs dz/(l + Zobs) then we have 

Ni s o- lsc (v) dv 



1"c(Ze 



. Vobs) = C I 



H 



(3) 



Fixing the emission frequency and the emission redshift one can 
immediately determine the current redshift as a function of v 
with v/v em = (1 +z)/(l +Zem)- Now the escape probability in the 
Lyman-continuum can be directly given in terms of the normal- 
ized emission profile = No-\ sc {v) Y 1 e - hv l k ^ T y M 



^Ly-c 



fern) 



-r 



,rec 
■em 



(v') exp [-T c (z e 



c)] dv', (4) 



where N is determined from the condition J (/^(v) dv = 1. 
Because 0^ is sufficiently narrow (Av/v c ~ k^Tylhv c ~ 1 /40 at 
z ~ 1400) it only takes a very short time to cross the emission 
profile via redshifting. Since the characteristic time for changes 
of N\ s and H is much longer, one can write 



poo /-•) 
<Jv c |_ kJv c 



c(v) 
dv 



dv' (5) 



where fruciv) = cr lsc (v)/cri sc (v c ) and the optical depth at the 
threshold frequency is To(z em ) = c o- lsc (v c ) Ni s (z em )/ H^). 

It is possible to simplify Eq. © even further. First, since 
the emission profile is very narrow one may use <f>f^ as 
ft e -h[v-v c yk s T y ^ w j[j 1 ]- ne corresponding normalization constant 
N. Assuming <xi sc m [v c /v] 3 one then obtains 

1 



3 Ly-c 



(Zem) 



f 
Jo 



-eii+if 



1 + T^ SC ' 



(6) 



where r^ sc = To^p 6 . We see that interpreting the continuum 

profiles as a very narrow lin^E P^Ic ° appears to be similar to 
the Sobolev approximation for the escape of photons from res- 
onance lines. In particular P^ c c also scales as 1/t^ sc for large 
optical depth r[r sc » 1 . However, unlike in the Sobolev approx- 
imation there is no symmetry between the absorption and emis- 
sion profile, which actually is one of the key requirements that 
makes the Sobolev escape probability completely independent 
of the line profile. Comparing the results of Eq. (0 and Eq. © 
shows that the latter equation provides a sufficient approxima- 
tion to p\l c c (see Fig.|5]l. 

3 Here and below we neglect the difference in the electron and photon 



temperature, which is smaller than AT/T ~ 10" 
during recombination. 
4 This is possible due to the exponential factor in ( 



10 4 at most times 
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Using Eq. © it is easy to include the possibility of direct 
recombinations to the ground state connected with the escape of 
continuum photons. For this, one should add 

Afl cIs * Pejrfem) X [N ls R 1s o ~ N e N p R cls ] , (7) 

to the rate equations of the electrons and subtract this term for the 
ls-equation. Here R c i s and R[ SC are the recombination and pho- 
toionization rate of the ground state, respectively. Note that we 
assume that on the blue side of the Lyman-continuum the spec- 
trum is Planckian. Therefore the photoionization rate R\ sc does 
not include any spectral distortion. However, due to the release 
of photons during He I recombination some small distortions are 
created, which lead to an increase of the photoionization rate. 

Comparing the escape probability for the Lyman-continuum 
with the Sobolev escape probability for the Lyman-*? photons 
(see Fig. |5]l one can see that p\l c c (z e m) is roughly 10 - 100 times 
larger at most times. However, including the possibility of direct 
recombinations to the ground state in our recombination code 
yields a tiny correction to X e (see Fig.|6]l, which shows that AR C [ S 
as given by Eq. (0 is still many orders of magnitude smaller than 
the net Lyman-o- transition rate. One can safely neglect direct 
recombinations to the ground state of hydrogen for computations 
of the recombination history. 

4. Discussion 

We have shown that Lyn — > Ly(n - 1) feedback leads to a modifi- 
cation in the ionization history of AN e /N e ~ 0.22% at z ~ 1050. 
Since this is very close to the peak o f the Thomson-visibility 
function (ISunvaev & Zeldovichll 19701) . corrections to the CMB 
power spectra of similar order are expected. Computing the 
CMB TT and EE power spectra using the results for the ion- 
ization history with and without the additional feedback we find 
changes with a typical peak to peak amplitude \ACf /Cf \ ~ 
0.15% and |ACf £ /Cf £ | ~ 0.36% at l < 3000 (see Fig.Q. 

As shown by Chluba & Sunvaevl d2007l) due to two-photon 
processes one can expect some asymmetries in the emission pro- 
files of the higher Lyman-series lines. This could modify the 
amount and time-dependence for the Lyn — > Ly(n - 1) feedback 
process and may lead to additional differences in the results pre- 
sented above. However, because the bulk of photons that escape 
in the red wing of a particular Lyman-transition is expected to 
come from close to the line center, this may be of minor impor- 
tance. Here the most interesting aspect may be the changes in 
the time-dependence of the feedback, but more detailed compu- 
tations would be necessary to understand this problem. 

In addition one should look at the feedback-induced cor- 
rections in the hydrogen recombination history due to the He I 
spectral distortions, in particular due to the photons appearing 
in the 2'Pi — > 1 'So resonance transition, the 2 3 Pi — > l'Soin- 
tercombination line and the 2 'So — » 1 'So two-photon contin- 
uum. It is clear that all high frequency He n lines will be re- 
processed during He I recombination, and correspondingly can 
only affect the ionization history during that epoch. However, 
for the feedback of the He I lines on hydrogen recombination it 
will be important to compute the re-processing of photons in the 
Lyman-continuum and all the subsequent Lyman-series transi- 
tions. Since all of these transitions are very optically thick a huge 
part of the photons released during He I and He II recombination 
will never reach the observer today, but due to re-absorption by 
neutral hydrogen, they are fully convert ed to hydrogen Lyman-q 
and 2s- Is continuum photons (see a lso Kholup enko et al. 2007; 
Rubino-Martm et al. in preparation). Since most of the helium 



2'Pi l'So photons are released at z ~ 1800 - 2600, they 
should be re-absorbed by neutral hydrogen atoms at the early 
stages z > 1200 - 1400. One can also conclude this from 
the paper of Kholupenk o et alj (|2007), where the re-processed 
He II — > He I photons appear on the red side of the hydrogen 
Lyman-ff distortion (most hydrogen Lyman-c photons are re- 
leased at z ~ 1400). Therefore the helium photons affect the hy- 
drogen recombination history well before the maximum of the 
Thomson-visibility function (z ~ 1100), and hence should have 
a rather small impact on the CMB power spectra. However, a 
more careful comp utation is required and will be desc ribed in a 
forthcoming paper (Rubino-M artm et alJIin pr eparation). 

Acknowledgements. The authors thank the anonymous referee for his useful 
comments. 
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Fig. 4. Changes in the net 2s- Is two-photon decay rate due to the 
inclusion of Ly« — > Ly(« - 1) feedback, as labeled. The results 
are based on the computation for a 15-shell hydrogen atom. 
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Fig. 5. Escape probability for the Lyman-continuum as given by 
Eq. (0 and Eq. ([§}. For comparison we also show the Sobolev 
escape probability for the Lyman-a- photons. 
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Fig. 3. Feedback-induced changes of A/ V (z e m) for the Lyman-a- 
and the modifications of the Lyman-/? and y lines before their 
absorption within the corresponding lower lying resonance. The 
results are based on the computations for a 15-shell hydrogen 
atom. 
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Fig. 6. Changes in the free electron fraction due to inclusion of 
Lyman-continuum escape relative to the reference model with- 
out any direct recombinations to Is. The computations were per- 
formed for a 5-shell atom where for the Lyman-series escape the 
Sobolev approximation was used. Note that the absolute value 
of AXJX e is ~ 1(T 6 only. 




Fig. 7. Changes in the CMB temperature (TT) and polarization 
(EE) power spectra. The diffe rences were computed using our 
modified versions of Cmbeasy (Doran 2005), which allows load- 
ing of pre-calculated recombination histories, and where the cor- 
responding RECFAST-routine was improved to achieve higher nu- 
merical accuracy with solvers from the NAG-library. 



